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Abstract 



We develop an advanced method of solving homogeneous and inhomogeneous Bethe- 
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Salpeter equations by using the expansion over the complete set of 4-dimensional 
spherical harmonics. We solve Bethe-Salpeter equations for bound and scattering 
states of scalar and spinor particles for the case of one meson exchange kernels. 
Phase shifts calculated for the scalar model are in agreement with the previously 
published results. We discuss possible manifestations of separability for one meson 
exchange interaction kernels. 
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*sO : 1 Introduction and preliminaries 

Obtaining solutions of inhomogeneous Bethe-Salpeter (BS) equations is a long standing 
problem with clear motivation - BS amplitude in continuum (or half-off-shell T-matrix) 
is a necessary ingredient for description of final state interaction in numerous processes 
with the correlated nucleon-nucleon pair in final state. As an example one can mention, 
for instance, charge-exchange reaction on deuteron pD — > n(pp) pQ, deuteron break- 
up with fast forward pp-pair [2], electro- and photo-disintegration of deuteron [Hj etc. 
One of the possible ways to obtain the BS amplitude in continuum is the so-called one 
iteration approximation scheme [1], when the relativistic P-components are obtained as 
first iteration from the non relativistic wave function with a certain interaction kernel. 
In jj] the inhomogeneous BS equation for spinor particles is regularly treated with a 2- 
dimensional Gaussian mesh. Although there are no doubts in correctness of these results, 
it can be shown that, in general, inappropriate choice of a 2-dimensional mapping in 
integration procedure can lead to under- or overestimation of the solutions. In this sense 
a 2-dimensional change of variables that one has to apply in partial BS equation for 
transforming the limits of integration to finite ones is not well defined procedure. Thus, 
the search of the rigorous method to solve such an equations should be continued. 

In jH] the expansion of the BS amplitude and interaction kernel over hyperspherical 
basis is employed to formulate a method of solving the BS equation for bound states of 
scalar particles, 

m=ij -0^V(k,p)S( Pl )S( P2 )T(p), (1) 

Pi t 2 = P/2 ± p, and S(p) = (p 2 — m 2 + is)^ 1 . The corresponding expansions after the 
Wick rotation looks like 

T(ik 4 ,k) = $>r$)^nlmK), (2) 
nlm 
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where by definition k = yk 2 + k 2 , and = (%, 9, <fi) - angles of vector k = {k^k) in 
4-dimensional Euclidean space. Hyperspherical harmonics 

Z„i m (x, 0, 4>) = X nl (x)Y lm (0, 4>), X nl ( X ) ~ sin' X C l n t)(cos X ) , 

are proportional to the product of usual spherical harmonics for orbital momentum / and 
Gegenbauer polynomials C^t*. 

Performing a Wick rotation in ((TJ) and substituting © and (j3J) one reduces the BS 
equation to the infinite set of 1-dimensional integral equations. Thus, there is only one 
degree of freedom left in the choice of the mapping for integration. It is shown in [B], 
that this approach is very effective for numerical purpose. In particular, due to rapid 
convergence of hyperspherical expansion, the consideration of only first few equations 
from the set is sufficient. 



2 Inhomogeneous BS equation for scalar particles 

In this paper the method introduced in jB] is applied to the scalar inhomogeneous BS 
equation, and generalized to spinor equations, both homogeneous and inhomogeneous. 
Let us first consider the BS equation for the scattering of scalar particles: 

T(k) = V(k,q)+i j -0- 4 V(k,p)S( Pl )S(p 2 )T(p). (4) 

Here T(k) is half-off-shell T-matrix for the scattering of real particles with momenta 512, 
and in c.m.s., where the real 3-momentum of particles \q\ = p, 

s = (q l + q 2 ) 2 = (2E i> ) 2 , E p = ^m 2 +p 2 , q = (2m) 2 < s < (2m + /i) 2 

Before Gegenbauer decomposition of type Q and (j3J) the Wick rotation should be applied. 
It is more convenient to perform firstly in (JH) the usual partial expansion over the spherical 
harmonics Yi m (see also [Zj): 

Ti(k ,k) = g 2 Qi(k ,k;0,p) -i J "T^rf g 2 Qi{h,k- : p^p) S{p 1 )S{p 2 )T l (p ,p). (5) 

Here Qi(k ,k;pQ,p) = Qi(\k 2 + p 2 + fi 2 — (k — p ) 2 — ie]/2kp) is the adjoint Legendre 
function of the 2nd kind. The normalization of T-matrix is fixed by the free term in (JHJ) 
and leads to the following expression for the phase shifts: 

T z (0,p) = IQtx p^se iSl sin 5 h (6) 
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To eliminate the removable singularity at (po,p) = (0,p), we present the T-matrix in the 
following factorized form jH]: 



T^k) ~ tpiik^k)^), t,(s) ~ e iSl sin 5. 



thus obtaining from (jSJ) the equation for 

<pi{k ,k) = g 2 Qi(k ,k;0,p) -i J {g 2 Qi{k , k;p ,p) 

Qi(k Q ,k;0,p) Qi(0,p;p ,p)} S(p 1 )S(p 2 ) fi(po,p), 



Qi(s, 



(7) 



Qi(s) = Qi(0,p;0,p). In the following only the cases / = 0, 1 will be considered. After 
the Wick rotation the Gegenbauer decomposition of tpi for I = takes the form 

oo 

ip (ih,k) = kpJ2 9j( k ) X 2j ~ 2 ,o(x), (8) 

3=1 

and for the coefficient functions gj the system of integral equations is obtained: 
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The additional unknown function ro comes after the Wick rotation from the residue 
at the pole p 2 = \fs/2 — E p + ie, which is able to cross the imaginary p axis. It is 
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proportional to the T-matrix in Minkowsky space for certain value of relative energy, 
T o(p) ~ Tq{*/s/2 — E p ,p). The partial kernels are 

X m { X )X kl { X ) 



SlkiP) = / dxsm 2 X j^ 2 . .> . .-> , ■ 

j [p z — p z ) z -\- sp* cos 1 x 

77 1 
N ml (k,p) = J dxpsin 2 x P ^r{Qm(k^k;ip 4 ,p) 



+ Q m (-k 2 ,k;ip A ,p)} — 



2kp 

Xl m (Xp) 



(p 2 — p 2 ) 2 + sp 2 COS 2 Xp ' 

^ 1 

W n i(k,p) = J dxksm 2 XkX n i(xk) i~{Qi(ih,k;pl,p) + Qi(ik A , k;-p° 2 ,p)}, 
o ^ 

Ui{Kp) = ^{Qi(klk;p° 2 ,p) + Q l (k 2 ,k;-plp)}. 
It is easy to see that all the partial kernels are real expressions. 



3 Numerical results and separability 

As all the integrations in system ®, (|1U|) are 1-dimensional, the integrals can be replaced 
by finite sums by using the gaussian mesh. In this way this system is reduced to the usual 
system of linear equations, which can be solved by any appropriate method. Besides, 
the solution of integral equations ©, (fTUj) can be obtained by iterating the free term, 
thus constructing the Neumann series, which represents the solution of an inhomogeneous 
integral equation. It was explicitly checked that these two procedures of finding a solution 
lead to the same results. 

Taking the first five terms in the expansion (jHJ), the system Q, (fTUj) was numerically 
considered for the following set of parameters: 

■L- = 4vr, m = ii = 1 GeV, p = 0.77 GeV/c 
An 

In Fig. ^the numerical results for the first three coefficient functions from (jHJ) are shown. 
It was established that convergence of this expansion is quite rapid, that is why only 
3 components are presented and g 2 and g 3 are multiplied by 10 and 100 respectively. 
Centered symbols correspond to the points of the mesh, where the solution is defined. 
Lines connecting the points reproduce the results of the fit for the obtained Gegenbauer 
components with the fitting functions 



1. For gi{p),p =p 



* a]p 2 i- 2 



2. For g 2 (p) 



2 4 a 2 2j-2 

F(p) = -if—roV] , I , 2 ,. , (12) 
W P 2 + b 2 2 ^{p 2 + b 2 2 ^ V ; 
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3. For g 3 (p) 

.,2 I 2 4 „3„2j 

F(p) 



p 

p 2 + b 2 3 



y (13) 



12 3 

For the sake of brevity we don't show here the numerical values of the parameters ay ' 
and &i,2,3, which can be adjusted by any appropriate method (for instance, the built- 
in fitting procedure in Microcal Origin 7.0 can be employed). It should be mentioned 
that the quality of this analytical fit is excellent in the whole range of the argument p 
despite the very simple form of the functions This fact can be treated as an 

indication to some kind of separability in the one meson exchange interaction and will be 
discussed below. Besides, it was found that such a fit is valid for Gegenbauer components 
of solutions of any BS equation in ladder approximation (scalar and spinor equations, 
both homogeneous and inhomogeneous) . Hence, it is very useful for practical purposes, 
e.g. to represent the numerical results in a compact form. 

To discuss the manifestations of separability in more details, let us now turn to the 
homogeneous BS equation (JTJ). After Wick rotation and usual partial decomposition it 
gives 

/dp^ dp pp 1 
(27r) 3 Vi(ih,k;ip 4 ,p) S(p l )S(p 2 )Vi(ip A) p). (14) 

It is generally known that for the conventional separable kernel of the form (see e.g. [Oj) 

Vi(ik 4l k;ip 4 ,p) = J2 9i(k 2 ) g0) 

the vertex functions are expressed in terms of functions g: 

Ti(ik 4 , k) = J2 c i9i(k 2 ). 

i 

In our case, solving equation (fT4*|) for kernel (jSJ) by the method given above and fitting the 
calculated Gegenbauer components with analytical expressions like (|TT|) - (fT3^) . we obtain 
the following approximate representation for the solution (cf. (JSJ)): 

JV 

T^k) = Y J c)R J i(k 2 )X ]l { Xk ), (15) 

3=1 

where the fitting functions are 

R (~k 2 ) = ( * 2) "* 

It is clear that constructing the expression for the separable kernel like 

Vi(ik 4 ,k;ip A ,p) = 9u{p*,p) 9ji(h, k) (16) 

ij 

with 

g j i(k 4 ,k) = R jl (k?)X jl ( Xk ) (17) 
one identically reproduces the analytical form of the solution |T5|) . Therefore, the kernel 



pfijl . (fTTjl can be referred to as some general form for constructing separable kernels. Note 
that it has the following properties: 



5 



A dependence on mass (or s) in the functions g 

An explicit dependence on fc 4 in g coming from Gegenbauer polynomials being the 
functions of cos^fe = k±/k 

C 

At low k 2 (|16j) tends to a Yamaguchi kernel, because g(k 2 ) 



k 2 + p 2 



A concrete example of the solution of eq. (jT4*j) is given in Table 1. The parameters \j 
[GeV 5 ] for the partial kernel (jTUjl of rank 3 are 

An = 0.19613 A 22 = 0.13734 A 33 = -0.01644 
A 12 = -0.09675 Ai 3 = 0.00236 A 23 = -0.00034. 



J 




&o[GeV] 




V J0 





0.17904 


0.39736 
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2 


-0.00923 


0.41533 


1 


2 


4 


0.00127 


0.43362 


2 


3 


6 


-0.00055 


0.45127 


3 


4 


8 


0.00017 


0.47019 


4 


5 



Table 1. Numerical values of parameters in (J15|) for I = 0, m = 1 GeV, mass 
of bound state M = 1.9 GeV, /x = 0.1 GeV. Coefficients c° for odd values of j 
are equal to zero. 

Analyzing the obtained solutions for the system (JUJ), (jTUjl it is also informative to calculate 
corresponding phase shifts. The results for I = and I = 1 are presented in Figs. 121 and 
HU respectively, for different values of coupling constant A connected with g 2 as 

£ = 4vrA. 

4:71 

It is obvious, that for A = 0.7 there are no bound states in this system due to the 
Levinson's theorem. For the values A = 1, 3, 5 there is one bound state, and two bound 
states for A = 7. For I — 1 no bound states found for these values of A, as it follows from 
Fig. El Obtained results for the phase shifts are in a good agreement with jH]. 



4 Spinor BS equations 

Above the procedure of solving the BS equation and handling with its solutions is de- 
scribed for the case of scalar particles. Spinor equations are much more complicated, and 
so we do not show here the detailed formulas like (J0J), (jTUjl . Nevertheless, it should be 
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stressed that all the basic steps remain the same. Below we briefly mention only the main 
results for spinor case. 

The problem of obtaining the BS amplitude in the continuum (or half-off-shell T- 
matrix) for spinor particles as a solution of the inhomogeneous BS equation with realistic 
interaction kernel seems to be of great importance. We present the results of solving the 
BS equation in the 1 So channel 

f(k)=V(k, q )T(l)Tt s + (q)T(2)+i J ^ V(k,p) S( Pl ) f(p) S(p 2 ) T(2) (18) 

in the ladder approximation for the cases of scalar (r(l, 2) = 1) and pseudoscalar (r(l, 2) = 
75) one meson exchanges without cut-off. Here 

~ p + m ~ p-m 

S[P) = — 9 ^— , S(p) 



pi _ m 2 _|_ fa p2 _ m 2 _|_ £ £ 

In f!18|) the BS vertex function T(k) being the matrix 4x4 before Gegenbauer decom- 
position of type (J2J) and (jHJ) should be expanded over some full set of matrices. In 1 5'o 
channel there are only 4 such a matrices usually called spin-angular harmonics (for their 
explicit form see [HI]). 

r ls ++(-fc), r ls --(-fc), r 3 pe(-fe), r 3p? (-k), (19) 

and corresponding expansion after the Wick rotation reads like 

f{ik il k) = ^2t a (ik i ,\k\)r a (-k). (20) 

a 

Unfortunately, after the expansion (|2T)j) the Gegenbauer decomposition of resulting partial 
equations for t a (ik4, \ k\) can not be performed in a simple analytical way. To make it 
possible, another set of spin-angular matrices for 1 S'o channel was found: 

<T(U\ ^ 5 T(U\ T(U\ <T(U\ C01 \ 



and, similarly to (|2~Uj) . 



f(ik 4 , k) = ^2 9j(ik4, \k\)Tj(k). 
3=1 



It can be easily shown that there is a non-degenerate transformation from set (|2ip to 
(HHJ). The point of using these two bases in the 1 5' subspace is that the set (|2*T|) allows 
us a very simple Gegenbauer decomposition of partial equations for gj{ik^ \k\) in an 
analytical form, while the set is preferable in calculating the matrix elements of the 
real processes jTJ 0] , because of its clear physical meaning. 

From the statement above follows that the half-off-shell T-matrix T{k) is expanded 
over the basis (|2*T|) and the derived set of integral equations for the scalar coefficients 
after the Wick rotation is decomposed into Gegenbauer polynomials like in (}8|)- fT0|) . We 
solve this set of equations numerically, for m = 0.938 GeV, /x = 0.14 GeV, y/s = 2m + 
0.1 GeV, and here the final results for solutions are presented in the form of partial 
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components t a {ik±, \k\) of (J2Uj). In Figs. EHZ1 these functions are shown for the case of 
a scalar interaction, and in Figs. ISHTD for a pseudoscalar one. For scalar exchange the 
ratio between different components looks very familiar: t ls ++ is the largest one, the other 
components are much smaller. Note that this situation does not hold for pseudoscalar 
exchange: here ti s — becomes the main component, and the P-wave Upe is larger than 
tig++. Although such a ratio between components seems to be exotic, it is valid also 
for the bound state problem in the 1 So channel (eq. (|2*2|) with r(l,2) = 75). Moreover, 
it qualitatively reproduces the previous results (huge P-wave) for the BS amplitude in 
continuum obtained in [2 within the one iteration approximation scheme. 
Finally, the homogeneous spinor BS equation of the form 

*(*) = */ ^V(k,p)T(l)S( Pl )np)S( P 2)T(2), (22) 

was examined for one meson exchange interaction kernels with T(i),i = 1,2 as scalar, 
pseudoscalar, vector, tensor vertices. Only 1 S bound states were considered. In this 
study we followed the motivations of [TT], and the obtained results will be a subject of a 
separate publication. 

5 Conclusion 

The method described before was successfully applied to solve both, homogeneous and 
non homogeneous, BS equations for scalar and spinor particles. It shows a high efficiency 
and accuracy, and, in particular, does not require any extraordinary computer facilities. 
The obtained results are in a good agreement with calculations of other groups jS] and 
also agree with our previous results j2]. Also the manifestations of separability for one 
meson exchange kernel are discussed. From the analysis of the obtained exact solutions, 
a possible general form of the separable kernel is suggested. 
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|p|, GeV/c 

Figure 1: Gegenbauer components from (jHJ). Functions g% and are multiplied by 10 and 
100 respectively in view of their smallness. 
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Figure 2: Phase shifts (in the units of n) for I = and different values of coupling constant. 
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Figure 3: The same as in Fig. El for 1 = 1. 
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Figure 4: 1 Sq + partial component (the largest one) of the solution of eq. (jl8j) for scalar 
meson exchange interaction. 
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Figure 6: 3 P 1 6 partial component for scalar meson exchangt 
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Figure 7: 3 P° partial component for scalar meson exchange. 
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Figure 8: 1 Sq + partial component of the solution of eq. (jl8J) for pseudoscalar meson 
exchange interaction. 
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Figure 10: 3 P 1 e partial component for pseudoscalar meson exchange. 
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Figure 11: 3 P° partial component for pseudoscalar meson exchange. 
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